log using "$Rep_smokelabor/2_analysis/output_logs/appendix_figure2.log", replace
	
	use "$Rep_smokelabor/1_build/regdata/county_day.dta", clear
	
	tsset countyfip rfrnc_dt
	
	* Wind dir. FEs
	reghdfe pm25 F(20(-1)1).hms_deep hms_deep L(1(1)20).hms_deep, a(county_doy st_yr g60_wdir) cluster(COUNTY10 rfrnc_dt) 
	
	* State x wind dir. FEs
	reghdfe pm25 F(20(-1)1).hms_deep hms_deep L(1(1)20).hms_deep, a(county_doy st_yr fe_wdir_st) cluster(COUNTY10 rfrnc_dt) 
	
	* County x wid dir. FEs
	reghdfe pm25 F(20(-1)1).hms_deep hms_deep L(1(1)20).hms_deep, a(county_doy st_yr fe_wdir_county) cluster(COUNTY10 rfrnc_dt) 
	
	** plot
	import excel using "$Rep_smokelabor/2_analysis/output_figures/appendix_figure2.xlsx", clear sheet("pm25_winds") first
	
	tw 	connected b all_wind  eday, col(blue "237 125 49") msize(0.8 1) msymbol(circle X) || ///
		line st_wind county_wind eday, col(black teal%80) lp(dash solid)  ///
		xtitle("Days since smoke exposure") ytitle("Change in PM2.5 (ug/m3)") ///
		xline(0, lcol(black) lp(dot) lw(0.6)) ///
		yline(0, lcol(black) lp(dot) lw(0.6)) ///
		ylab(, nogrid) ///
		legend(region(lcol(white)) order(1 "Baseline specification" 2 "Wind dir. FEs" 3 "State x wind dir. FEs" 4 "County x wind dir. FEs")) ///
		graphregion(col(white))
	gr export "$Rep_smokelabor/2_analysis/output_figures/appendix_figure2.pdf", replace
	
log close 